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1.0  SUMMARY 


An  experimental  and  computational  investigation  of  supersonic  combustion  flows 
has  been  carried  out.  The  principal  objective  of  the  research  was  to  gain  a  more 
fundamental  understanding  of  mixing  and  chemical  reaction  in  supersonic  flows.  The 
research  effort  comprised  three  inter-related  elements:  (1)  stability  analyses  and 
numerical  simulations  of  compressible  reacting  flows;  (2)  development  of  laser-induced 
fluorescence  techniques  for  time-resolved  multidimensional  imaging  of  species 
concentration,  temperature,  velocity  and  pressure;  and  (3)  an  experimental  study  of 
mixing  and  combustion  in  a  supersonic  plane  mixing  layer  with  the  additional 
development  of  simple  mixing  enhancements.  The  specific  objectives  and  results  of  the 
research  of  each  of  these  program  elements  are  summarized  in  this  report. 

2.0  INTRODUCTION 

Air-breathing  propulsion  systems  offer  the  potential  of  higher  performance  than 
conventional  rocket  engines  for  hypersonic  flight.  To  realize  this  potential,  new 
combustor  design  concepts  are  required.  In  particular,  in  order  to  minimize  losses 
associated  with  strong  shock  waves  and  high  combustor  inlet  temperatures,  it  is  desirable 
to  maintain  high  flow  velocities  in  the  combustion  chamber.  This  design  concept  leads  to 
a  new  class  of  propulsion  devices  where  combustion  takes  place  in  supersonic  flow. 

Combustion  in  supersonic  flow  is  fundamentally  different  from  combustion  in  the 
subsonic  flow  regime  employed  in  all  currently-operating  aircraft  engines.  Many  of  the 
design  approaches  developed  over  the  years  for  subsonic  combustors,  e.g.,  ignition  and 
flame  stabilization  techniques,  are  not  applicable  to  supersonic  combustion  devices,  and 
the  current  understanding  of  the  fundamental  aspects  of  supersonic  combustion  is 
inadequate  to  support  development  of  these  devices. 

Recent  advances  in  diagnostic  capabilities  and  significant  improvements  in  our 
ability  to  compute  such  flows  offer  new  opportunities  to  obtain  needed  fundamental 
understanding  of  compressible  turbulent  reacting  flows.  To  achieve  this  understanding, 
we  have  carried  out  a  closely-coordinated  experimental  and  computational  program  that 
utilized  state-of-the-art  experimental  techniques  and  computational  methods. 

The  principal  objective  of  the  research  is  to  gain  a  more  fundamental 
understanding  of  the  flow  physics  and  chemistry  interactions  in  compressible  turbulent 
reacting  flows.  The  project  comprised  three  inter-related  efforts;  (1)  stability  analysis 
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and  numerical  simulation  of  supersonic  flows  with  mixing  and  chemical  reaction;  (2)  an 
experimental  study  of  mixing  and  combustion  in  supersonic  flows  with  the  additional 
development  of  simple  mixing  enhancements;  and  (3)  development  of  laser-induced 
fluorescence  techniques  for  time-resolved  multidimensional  imaging  of  species 
concentration,  temperature,  velocity  and  pressure  in  supersonic  flows.  A  close  coupling 
of  these  efforts  was  maintained  in  order  to  maximize  our  understanding  of  compressible 
turbulent  reacting  flows  with  an  emphasis  on  supersonic  combustion.  The  specific 
objectives  and  results  of  the  research  of  each  of  the  program  elements  are  described  in  the 
following  sections. 

3.0  STABILITY  ANALYSIS  OF  FLOW  STRUCTURE 

3.1  Objective 

The  objective  of  the  numerical  portion  of  the  project  has  been  to  use  stability 
analysis  to  gain  insight  into  how  the  flow  structure  and  physics  of  the  compressible 
reacting  mixing  layer  change  within  a  parameter  space  defined  by  compressibility,  heat 
release  and  the  density,  equivalence  and  velocity  ratios. 

3.2  Research  Results 

Previous  simulations  done  by  Planche  and  Reynolds  (1992a,b)  have  demonstrated 
that  a  compressible  reacting  mixing  layer  can  develop  two  instability  modes  in  addition 
to  the  more  common  central  mode  that  exists  unaccompanied  in  incompressible  non¬ 
reacting  flows.  These  two  additional  modes  are  termed  "outer"  because  of  their 
association  with  the  fast  and  slow  free  streams.  These  initial  simulations  demonstrated 
the  adverse  impact  of  an  outer  mode  structure  on  the  global  reaction  rate,  and  thereby 
motivated  our  interest  to  gain  better  insight  into  the  range  of  flow  conditions  that  result  in 
the  outer  modes  becoming  the  dominant  instability. 

3.2. 1  Explanation  and  Prediction  of  Flow  Structure 

The  density  weighted  vorticity  profile  CpdU/dy)  provides  an  extremely  useful 
tool  for  both  describing  and  predicting  the  first  order  behavior  of  the  mixing  layer.  A 
view  of  how  this  profile  changes  under  the  influence  of  increasing  compressibility  and 
heat  release  is  shown  in  Figure  1 .  Our  results  parameterize  the  effect  of  compressibility 
with  the  convective  Mach  number.  Me,  and  the  amount  of  heat  release  by  a  non- 
dimensional  heat  of  reaction,  0.  This  figure  demonstrates  that  both  combustion  (through 
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heat  release)  and  compressibility  (through  viscous  heating)  act  to  increase  the 
temperature  and  correspondingly  decrease  the  density  along  the  centerline.  This  results  in 
a  fundamental  change  in  the  ^  duldy  profile  from  a  single  centrally  located  peak  to  a 
bimodal  appearance  for  sufficiently  large  ©  and  Me  values.  This  evolution  suggests  that 
the  mixing  layer  will  develop  two  additional  instability  modes-one  associated  with  each 
of  the  fast  and  slow  streams-that  are  distinct  from  the  central  mode,  which  itself  is 
responsible  for  the  familiar  vortical  structure  found  in  incompressible,  non-reacting 
situations. 

The  implications  of  this  p  dU/dy  change  are  presented  in  Figures  2  and  3.  For  the 
case  of  low  compressibility  and  heat  release  shown  in  Figure  2,  the  single  peaked 
p  duldy  profile  results  in  a  flow  structure  that  provides  direct  entrainment  of  fuel  and 
oxidizer  together.  In  contrast  to  this,  when  both  compressibility  and  heat  release  are 
increased  in  Figure  3,  the  bimodal  density  weighted  vorticity  profile  results  in  a  two-step 
mixing  process:  the  vortical  structure  of  the  fast  mode  mixes  oxidizer  into  the  reaction 
zone  while  the  slow  mode  mixes  in  the  fuel.  Planche  and  Reynolds'  simulations 
documented  the  adverse  impact  of  this  two-step  process  on  the  global  reaction  rate, 
implying  a  reduction  in  the  mixing  performance  for  the  flow.  Given  this  impact,  there  is 
significant  interest  in  improving  our  understanding  of  how  the  structure  of  the  mixing 
layer,  and  hence  its  mixing  ability,  changes  as  a  function  of  the  full  parameter  space. 

3.2.2  Overview  of  the  Parameter  Space 

An  extensive,  five-dimensional  parametric  study  was  conducted  with  a  focus  on 
determining  whieh  variables  were  most  important  in  determining  the  flow  structure.  The 
results  indicated  that  compressibility,  heat  release  and  density  ratio  (s  effects 

were  the  most  important  parameters  in  establishing  the  flow  structure;  the  effect  of  the 
equivalence  and  velocity  ratios  were  found  to  be  less  critical  (Day  et  al.  1997b).  A  good 
appreciation  for  the  effect  of  the  three  most  important  parameters  can  be  developed  by 
considering  the  evolution  of  the  eigenfunction  energy  profile,  =  uu  +  vv*  +  ww* .  To 
create  a  visualization  that  combines  a  view  of  both  the  distribution  of  vortical  energy  and 
the  amplification  rate,  we  define  a  new  variable 

where  the  energy  term  is  calculated  for  the  instability  wave  at  the  frequency  and 
obliquity  angle  of  the  most  amplified  mode.  The  energy  was  normalized  to  have  a 
maximum  value  of  unity  before  being  multiplied  by  the  ratio  of  the  amplification  rate  at 
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the  position  of  interest  (-a,,max)  to  a  reference  amplification  (-a,,max,o)-  The  reference 
chosen  here  is  the  canonical  case  of  an  incompressible,  non-reacting,  equal  density 
mixing  layer  with  tJa  =  0 . 

The  'P  profile  that  corresponds  to  this  reference  condition  is  shown  in  Figure  4  as 
the  first  of  three  images  that  detail  the  influence  of  compressibility  on  the  flow  structure. 
There  is  only  a  single  amplified  mode  at  the  reference  condition  and  its  wide  and 
centrally  located  energy  profile  identifies  it  as  the  central  mode.  The  broad  distribution  of 
energy  in  this  case  is  the  signature  of  the  well  understood  large  scale  structure  that 
directly  mixes  fuel  and  oxidizer  together.  As  compressibility  increases  to  Me  =  1  in 
Figure  4,  the  growth  rate  of  the  central  mode  drops  and  its  energy  distribution  becomes 
slightly  more  compact.  Also  appearing  at  this  compressibility,  but  at  very  low 
amplification  rates,  are  the  fast  and  slow  outer  modes.  Consistent  with  the  ideas  presented 
in  our  p  dU/dy  discussion,  the  results  at  Me  =  2  indicate  that  the  growth  rates  of  the  outer 
modes  strengthen  relative  to  the  central  mode  with  increasing  compressibility,  but  still 
remain  secondary  to  the  rate  of  the  central  mode  at  this  Me  value. 

The  outer  modes  will  dominate  the  flow  structure  when  heat  release  effects 
become  significant.  This  situation  is  documented  in  Figure  5  which  shows  the  effect  of 
increasing  heat  release  relative  to  the  non-reacting  case  at  Me  =  2.  A  comparison  between 
the  eigenfunction  energy  profiles  at  Me  =  2,  0  =  0  and  Me  =  2,  0  =  1  demonstrates  that 
the  introduction  of  heat  release  at  this  Mach  number  entirely  stabilizes  the  central  niode 
but  destabilizes  the  outer  modes.  The  nature  of  the  two-step  mixing  process  that  was 
detailed  in  the  DNS  studies  of  Planche  and  Reynolds  (1992a,b)  is  evident  from  the 
separate  and  distinct  energy  distributions  shown  in  the  Me  =  2,  0  =1  plot. 

The  influence  of  the  density  ratio  on  the  flow  structure  is  also  seen  in  Figure  5. 
This  sequence  of  energy  profiles  details  the  effect  of  a  fast  (s  =  0.5)  and  slow  (s  =  2.0) 
side  density  bias  relative  to  the  equal  density  Me  =  2,  0  =  0  case.  Doubling  the  density  on 
the  fast  side  causes  the  fast  and  central  modes  to  be  nearly  equally  amplified,  forming  one 
variety  of  colayer  structure  in  the  flow,  while  the  same  change  is  observed  to  stabilize  the 
slow  mode.  The  increase  in  density  on  the  fast  side  has  also  moved  the  position  of  all 
three  vortical  modes  upwards  in  the  layer.  The  profiles  for  a  slow  side  density  bias  show 
the  opposite  result  with  this  change  causing  the  formation  of  a  central/slow  mode  colayer 
and  a  general  shift  towards  the  slow  side. 
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Fig.  1 :  Density  weighted  vorticity  development  with  compressibility  and  heat  release  (1/2  =  0,  s  =  1,  (/>  =  1). 


(a) 

Fig.  2:  Characteristics  of  a  low  compressibility,  low  heat  release  mixing  layer  at  0  =  0.1  and  Me  =  0.1:  (a) 
density  weighted  vorticity  profile  (b)  single  step  mixing  process  for  this  central  mode  flow  structure. 


(a) 


Fig.  3:  Characteristics  of  fully  compressible  reacting  mixing  layer  at  0  =  2,0  and  Me  =  1.0:  (a)  density 
weighted  vorticity  profile  (b)  two  step  mixing  process  for  this  outer  mode  flow  structure. 
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Fig.  6:  Evolution  of  the  modal  structure  regime 
chart  in  the  (s,0)  plane  with  compressibility 
through  cases  of  (a)  low,  (b)  transitional  and  (c) 
high  Me  values  ((/>  =  1, 1/2  =  0). 


Increasing 
Density  Ratio 
(s=P2/Pi) 


Fig.  7:  Evolution  of  the  modal  structure  regime 
chart  in  the  (Mc,0)  plane  with  density  ratio 
value  increasing  through  cases  of  (a)  s  =  0.5, 
(b)  s  =  1  and  (c)  s  =  2  (0  =  1,  U2  =  0). 
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3.2.3  Structural  Regime  Charts 

A  more  complete  understanding  for  the  interaction  of  the  three  dominant 
parameters  (Me,  ©,  and  s)  can  be  developed  by  taking  two-dimensional  cuts  through  the 
parameter  space  to  investigate  regions  where  specific  modes  will  dominate.  A  series  of 
these  regime  charts  are  presented  in  Figures  6  and  7  for  cuts  in  the  (s,  0)  and  (Me,  0) 
planes,  respectively. 

The  first  point  of  interest  concerns  the  shape  of  the  colayer  boundaries  for  low  Me 
values  in  (s,  0)  plan  of  Figure  6(a).The  evolution  of  the  p  du/dy  profile  toward  two 
outer  peaks  with  higher  heat  release  is  reflected  in  the  progression  of  the  regime  chart 
from  a  central  to  outer  mode  structure  with  increasing  0.  Also,  the  correlation  of  regions 
where  an  outer  mode  (fast  or  slow)  dominates  with  regions  where  their  companion 
streams  have  greater  density  follows  from  predictions  based  on  the  p  du/dy  profile.  The 
wedge  shape  of  the  central  mode  region  in  Figure  6(a)  can  be  explained  from  our  earlier 
observations  from  the  eigenfunction  energy  profile.  The  effect  of  increasing  heat  release 
to  stabilize  the  central  and  destabilize  the  outer  modes  is  reflected  in  the  diminishing 
density  ratio  range  where  the  central  mode  dominates  as  heat  release  is  increased;  the 
outer  modes  require  less  bias  in  the  density  ratio  value  to  match  the  central  mode 
amplification  as  the  0  value  rises.  The  conditions  of  particular  interest  in  this  regime 
chart  are  the  locus  of  points  that  form  the  boundaries  between  modes.  For  these 
conditions,  the  flow  will  develop  two  vortical  modes  that  contribute  to  the  mixing  of  the 
fast  and  slow  streams  and  therefore  hold  promise  for  enhanced  mixing  performance. 

Shifting  focus  now  to  the  results  at  Me  =  0.5  that  also  appear  in  Figure  6(a),  we 
see  that  the  predominant  features  of  the  previous  chart  remain  unchanged  but  here  the 
outer  modes  occupy  a  larger  region.  This  follows  directly  from  the  fact  that  both 
compressibility  and  heat  release  generate  conditions  where  outer  modes  dominate  the 
structure.  Further  increases  in  Me  are  pictured  in  Figures  6(b)  and  (c)  and  they  show  a 
continuation  of  the  trend  toward  outer  mode  domination  and  three-dimensionality  for  all 
modes.  Indeed,  at  Me  =  2  it  becomes  possible  to  maintain  a  colayer  structure  in  non¬ 
reacting  flow,  provided  the  density  ratio  is  carefully  chosen. 

Figure  7  offers  a  perspective  from  the  (Me,  0)  plane  for  a  variety  of  density  ratios. 
These  planar  cuts  allow  discussion  of  the  mixing  layer  development  from  the  reference 
condition  of  an  incompressible,  non-reacting  case  (located  in  the  lower  left  comer)  to  a 
fully  compressible  reacting  flow  over  a  range  of  fixed  density  ratios.  It  is  useful  to  start 
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the  discussion  with  the  reference  case  of  s  =  1  that  is  shown  in  Figure  7(b).  Under  the 
conditions  of  Me  =  0  and  ©  =  0  the  familiar,  centrally  located,  two-dimensional,  large 
scale  structure  will  be  present  and  dominant.  With  increasing  compressibility  there  is  first 
a  transition  to  three  dimensionality  and  then,  for  Mg  >  2.5,  the  central  mode  gives  way  to 
the  slow  three-dimensional  mode  as  the  dominant  flow  structure.  With  increasing  heat 
release  the  progression  from  central  to  slow  mode  structure  is  again  observed.  At  this 
unity  density  ratio  value  the  slow  mode  will  slightly  out-pace  the  fast  and  it  dominates 
the  regions  of  both  high  compressibility  and  heat  release. 

For  the  case  of  a  fast  side  density  bias  of  s  =  2/3  pictured  in  Figure  7(a),  the 
central/outer  mode  colayer  boundary  maintains  the  same  characteristics  as  in  the  s  =  1 
case  but  here  the  fast  mode  replaces  the  slow  as  the  dominant  outer  structure  over  most  of 
the  considered  regime.  The  s  =  1/2  case  also  shown  in  Figure  7(a)  demonstrates  that 
colayer  regions  are  almost  entirely  eliminated  at  this  level  of  density  bias.  Therefore, 
there  will  only  be  one  discernible  mixing  mode  in  the  flow  for  the  vast  majority  of  this 
parameter  plane.  When  the  density  bias  is  moved  to  the  slow  side  of  the  layer  in  Figure 
7(c),  the  location  of  the  central/slow  colayer  drops  in  both  Me  and  0  values  and  also 
reduces  in  length  compared  to  the  equivalent  density  case  of  Figure  7(b).  When  the 
density  ratio  is  increased  further  to  s=2,  the  colayer  boundary  is  all  but  eliminated. 

4.0  DIAGNOSTIC  TECHNIQUE  DEVELOPMENT 
4.1  Objectives 

Over  the  past  several  years,  investigators  in  our  laboratories  have  pursued  the 
development  of  laser-based  diagnostic  techniques  applicable  in  the  pulsed  flow  facilities 
frequently  used  in  ground-based  testing  of  high-speed  aerodynamic  and  propulsion 
concepts.  Our  efforts  have  focused  on  PLIF  imaging,  which  is  particularly  attractive 
because  it  can  yield  species-  and  quantum-state-specific  data  with  high  temporal  and 
spatial  resolution,  excellent  background  rejection,  and  low  noise.  The  use  of  PLIF  in 
transient,  high-speed  flows  was  pioneered  at  Stanford  and  workers  in  our  laboratories 
were  able  to  demonstrate  essentially  instantaneous  measurements  of  rotational  and 
vibrational  temperature  and  velocity  fields. 
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4.2  Research  Results 


4.2.1  Flow  Facilities 

The  small-scale  facility  shown  schematically  in  Fig.  8  was  constructed  to  provide 
inexpensively  a  controlled  environment  for  PLIF  diagnostic  development.  The  facility 
may  be  operated  either  as  a  shock  tube  or  as  a  reflection-type  shock  tunnel  (with  a  rapid 
conversion  procedure)  and  may  be  used  to  create  a  variety  of  supersonic  flowfields.  A 
wide  range  of  operating  conditions  is  available  for  tests  in  the  facility.  Figure  9  illustrates 
the  shocked-gas  pressure  and  temperature  in  the  driven  section  as  a  function  of  initial  test 
gas  pressure  for  a  typical  driver-gas  condition.  Usually,  the  test  gas  is  composed  of  an 
inert  gas  {i.e.,  Ar  or  N2)  seeded  with  a  small  mole  fraction  of  a  fluorescent  species  (NO) 
or  with  reactants  that  form  a  fluorescent  species  (OH)  following  shock-heating. 
Alternatively,  NO-  or  reactant-bearing  gas  may  be  introduced  through  a  wall-mounted 
injector.  The  duration  of  steady  flow  conditions  varies  inversely  with  the  enthalpy  of  the 
flow  and  ranges  from  ~150  ps  to  1.5  ms. 

The  basic  optical  arrangement  used  in  acquiring  PLIF  images  is  illustrated  in  Fig. 
8.  Frequency-doubled  photons  from  a  tunable,  excimer-pumped  dye  laser  are  formed 
into  a  sheet  ~50  mm  wide  and  0.2-0.4  mm  thick  and  directed  across  the  midplane  of  the 
test  section  through  a  UV-grade  fused  silica  window.  During  each  run  of  the  facility,  a 
single  laser  pulse,  fired  at  the  desired  moment  of  observation,  excites  an  isolated 
rovibronic  transition  of  the  tracer  species  present  in  the  flowfield.  A  fraction  of  the 
resulting  fluorescence  burst  is  imaged  onto  a  gated,  intensified  charge-coupled  device 
array  mounted  perpendicular  to  the  laser  sheet.  Spectral  filtering  is  used  to  eliminate 
detection  of  luminosity  from  the  flow  as  well  as  fluorescence  resonant  with  the  laser 
frequency  and  direct  laser  scattering  from  walls  and  windows.  The  captured  PLIF  image 
is  then  transferred  to  computer  for  further  analysis  and  display. 

4.2.2  PLIF  Imaging  of  Shock  Tunnel  Transient 

Figure  10  shows  the  results  of  a  series  of  experiments  in  which  PLEF  imaging  of  NO  was 
used  to  examine  the  evolution  of  the  flow  up-  and  downstream  of  the  shock-tunnel  nozzle 
(Palmer  and  Hanson,  1995).  The  images  provide  midplane  views  of  the  reservoir  and  jet 
flowfields  at  different  delays  with  respect  to  the  time  at  which  the  incident  shock  reached 
the  endwall.  The  initial  distortion  of  the  reflected  shock  due  to  the  nozzle  entrance  in  the 
endwall  and  the  resultant  perturbations  in  the  reservoir  are  observable.  The  process  of  jet 
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startup  to  a  quasi-steady  flow  condition,  including  the  bursting  of  the  diaphragm  placed  at 
the  nozzle  exit  and  its  passage  through  the  downstream  imaged  region,  is  illustrated  as 
well.  At  long  delay  times,  disturbances  arriving  from  upstream  altered  the  pressure  and 
temperature  of  the  gas,  thereby  terminating  the  period  of  stable  reservoir  conditions. 
These  images  not  only  allow  a  detailed  examination  of  the  flow  physics,  but  also  provide 
crucial  information  about  the  operating  characteristics  of  the  shock  tunnel. 

4.2.3  PLIF  Imaging  of  Injector  Flowfield 

Because  the  PLIF  signal  is  a  complex  function  of  temperature,  pressure,  and  gas 
composition,  the  interpretation  of  an  individual  image  is  often  difficult.  However, 
analysis  of  two  or  more  PLIF  images  of  a  flowfield  can  yield  a  specific  flow  property, 
such  as  temperature  or  velocity.  For  instance,  investigators  in  our  laboratory  have  used 
two  complete  laser/camera  systems  in  the  acquisition  of  temporally-resolved  and  mean 
rotational  temperature  fields  in  a  model  SCRAMJET  injector  flowfield  (McMillin  et  ai, 
1994),  as  shown  in  Fig.  11.  The  lasers  were  tuned  to  two  transitions  originating  from 
different  rotational  states  within  the  lowest  vibrational  level  of  NO  or  OH.  A  brief  time 
period  between  the  firing  of  the  two  lasers  allowed  each  camera  to  integrate  the  signal 
induced  by  only  one  of  the  lasers,  while  capturing  the  flowfield  with  flow-stopping 
temporal  resolution.  The  rotational  temperature,  assumed  to  be  identical  to  the 
translational  temperature,  was  inferred  from  the  ratio  of  the  two  signals  on  a  pixel-by- 
pixel  basis.  A  comparison  of  the  results  obtained  with  PLIF  of  NO  and  OH  illustrates  the 
complementary  nature  of  the  two  techniques  and,  in  the  region  of  overlap,  validates  each 
of  the  mean  temperature  measurements. 

4.2.4  PLIF  Thermometry  in  Underexpanded  Free  Jet 

An  NO  PLIF  thermometry  technique  has  also  been  applied  in  a  highly  underexpanded 
free  jet  created  at  the  exit  of  a  converging  nozzle  (Palmer  et  ai,  1996)  .  The  results  of 
these  experiments  are  shown  in  Fig.  12.  Whereas  previous  rotational  temperature 
measurements  were  performed  in  flows  with  shock  waves,  the  focus  in  this  experiment 
was  on  application  in  a  high-speed  flow  with  yet  more  extreme  temperature  and  pressure 
variations.  The  method  used  was  similar  to  that  described  above,  although  only  one  PLIF 
image  was  acquired  during  each  run  of  the  facility.  Because  the  flow  was  laminar  and 
reproducible,  temporal  resolution  was  not  required,  expect  insofar  as  each  image  had  to 
be  captured  during  its  brief,  quasi-steady  period,  which  made  necessary  the  use  of  a 
single-shot  imaging  technique.  PLIF  images  acquired  by  pumping  transitions  originating 
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from  four  different  rotational  states  were  combined  to  yield  the  mean  temperature  field  in 
the  free  jet.  The  temperature  inferred  from  single-shot  and  frame-averaged  PLIF  images 
agrees  well  with  the  results  of  a  simulation  of  the  flowfield  by  the  method  of 
characteristics. 

The  examples  of  PLIF  results  provided  here  demonstrate  the  frame-on-demand 
capability  crucial  for  imaging  experiments  in  impulse  flow  facilities  as  well  as  some  of 
the  flow  property  measurement  strategies  developed  in  our  laboratory.  Other  applications 
of  PLDF  imaging  of  NO  in  the  shock  tube/tunnel  facility  have  included:  visualization  of 
the  establishment  of  bow  shocks  over  blunt  bodies  placed  in  the  shock  tube  (McMillin  et 
al,  1991)  and  at  the  exit  of  a  converging-diverging  shock-tunnel  nozzle  (Houwing  et  al, 
1996)  along  with  rotational  thermometry  in  the  resulting  quasi-steady  flows,  vibrational 
thermometry  across  planar  shock  waves  (McMillin  et  al,  1989)  and  in  free  jets  with 
vibrational  nonequilibrium  (Palmer  and  Hanson,  1993a),  and  temporally-resolved,  two- 
component  velocimetry  in  a  free  jet  flow  (Palmer  and  Hanson,  1993b).  In  addition,  PLIF 
imaging  of  OH  has  been  applied  to  measurements  in  a  variety  of  reacting,  high-speed 
flowfields,  including:  temporally-resolved  rotational  thermometry  in  planar  shocks 
(Seitzman  et  al,  1994)  and  in  bow  shock  flows  (Seitzman  and  Hanson,  1993), 
visualization  of  the  development  of  a  free  jet  of  combustion  product  gases  (Palmer  and 
Hanson,  1995),  and  mean  rotational  thermometry  in  the  resulting  quasi-steady  flowfield 
(Palmer  and  Hanson,  1996) . 
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Fig.  8.  Shock  tube/tunnel  facility  with  associated  optical  and  electronic  components  for  PLEF  imaging.  In 
the  configuration  shown,  the  facility  is  operated  as  a  reflection-type  shock  tunnel;  for  operation  as  a  shock 
tube,  the  nozzle  plate,  located  at  the  end  of  the  driven  section,  is  replaced  with  a  solid  endwall. 


Fig.  9.  Shock  tube  performance  characteristics  with  an  initial  driven-section  loading  of  (a)  N2  or  (b)  Ar  at 
pressure  Pjnit  and  300  K  for  a  single  driver-gas  condition  (/.^.,  room-temperature  He  at  3.4  atm,  with  a  10- 
mil-thick  mylar  diaphragm).  The  pressure  (P)  and  temperature  (T)  behind  the  incident  and  reflected  shocks 
were  calculated  using  incident  shock  speeds  measured  at  various 
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Fig.  10.  PLIF  imaging  of  the  shock-tunnel  flow  transient.  The  locations  of  the  two  imaged  regions 
reservior  images  and  jet  images)  are  illustrated  in  (a),  which  shows  the  test  section  midplane  viewed  from 
above  and  gives  the  nominal  conditions  in  the  reservoir  and  far  downstream  of  the  nozzle.  Each  NO  PLIF 
image  in  (b.1-20)  was  acquired  at  the  delay  times  listed  with  it  (with  respect  to  the  time  of  incident  shock 
reflection).  The  sequences  illustrate  the  process  of  shock  reflection,  the  development  of  the  free  jet 
flowfield,  and  the  establishment  and  duration  of  steady  flow  conditions. 
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(b)  NO  and  OH  PLIF  thermometry  results 
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Fig.  11.  PLIF  thermometry  in  a  model  SCRAMJET  injector  flowfield.  The  test  section  midplane,  viewed 
from  above,  is  shown  in  (a)  to  illustrate  the  locations  of  the  imaged  regions  (i.e.,  NO  images  and  OH 
images)  with  respect  to  the  step  in  the  side  wall  and  the  jet-injection  port.  The  images  in  (b.l)  are 
temporally-resolved  temperature  fields  acquired  with  PLIF  of  NO  and  OH,  as  labeled.  The  images  in  (b.2) 
are  mean  temperature  fields  constructed  by  frame- averaging  individual  images  of  the  type  shown  in  (b.l). 
Six  single-shot  images  were  used  for  NO  and  twelve  for  OH  in  (b.2).  Also  shown  is  a  composite  image 
constructed  from  the  mean  NO  and  OH  temperature  images.  The  scale  located  at  the  bottom  of  (b.l)  gives 
the  relationship  between  the  displayed  color  and  the  rotational  temperature  in  the  images. 
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(a)  test  section  midplane 


(b)  NO  PLIF  thermometry  results 
(1)  single-shot 


(2)  frame-averaged 
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Fig.  12.  Multi-line  PLIF  thermometry  in  a  shock-tunnel  free  jet  The  test  section  midplane,  viewed  from 
above,  is  shown  in  (a)  to  illustrate  the  location  of  the  imaged  region.  Shown  in  (b)  are  temperature  fields 
inferred  from  (1)  single-shot  and  (2)  frame-averaged  PLIF  images.  The  temperature  predicted  by  a  method- 
of-characteristics  simulation  of  the  free  jet  is  given  for  comparison  in  (b.3). 
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5.0  MIXING  AND  REACTION  IN  SUPERSONIC  FLOW 


5.1  Objectives 

The  objectives  of  this  part  of  the  program  are  (1)  the  experimental  study  of  mixing  and 
chemical  reaction  in  a  compressible  mixing  layer  and  (2)  investigation  of  mixing 
enhancement  in  supersonic  shear  layers.  The  data  acquired,  using  conventional  and  laser- 
based  diagnostics,  provide  important  new  insights  into  mixing  and  reaction  processes 
under  compressible  conditions. 

5.2  Research  Results 

In  the  primary  experiments,  the  convective  Mach  number  was  varied  from  Me  = 
0.25  to  0.75,  spanning  the  range  from  low  to  moderately  high  compressibility,  and  the 
influence  on  molecular-level  scalar  mixing  measured  with  the  relatively  new  cold 
chemistry  planar  laser-induced  fluorescence  technique.  The  technique  provided  a 
resolution-insensitive  measurement  of  pure  fluid  from  which  the  time-averaged 
probability  of  mixed  fluid,  shear  layer  mixing  efficiency,  mixed  fluid  composition,  and 
other  scalar  quantities  were  derived  from  flip  experiments  where  the  constituent  gases  of 
each  stream  are  exchanged.  More  traditional  passive  scalar  PLIF  measurements  were  also 
made  to  determine  the  mean  composition  profile  and  visualize  the  mixture  fraction  field. 

In  a  secondary  set  of  experiments,  the  use  of  boundary  layer  disturbances  on  the 
high-speed  splitter  tip  was  examined  as  a  means  to  enhance  mixing  with  low  pressure 
loss.  The  shape,  number,  obliquity,  offset,  and  thickness  of  the  disturbances  were 
parametrically  varied  and  the  qualitative  effect  on  mixed  fluid  visualized  with  planar 
laser  Mie  scattering.  For  the  most  promising  geometry,  quantitative  measurements  were 
made  using  the  cold  chemistry  mixing  diagnostic. 

Key  results  from  the  1)  diagnostic  development  effort,  2)  influence  of 
compressibility  on  scalar  mixing,  and  3)  mixing  enhancement  studies  are  summarized 
below.  The  complete  findings  of  this  work  are  described  in  the  Ph.D.  dissertation  of  T. 
Island  (1997). 
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Fig.  13.  Calculation  comparing  the  cold  chemistry  and  passive  scalar  signals.  The  seeded  high-speed  stream 
(^  =  1)  is  N2  and  the  unseededstream  is  air.  The  calculation  is  performed  with  P  =  10  psia,  =  150  K  and 
T2  =  290  K,  matching  the  high-compressibility  experimental  conditions  and  assumes  IM21  =  202  ns,  = 
0.01  and  =  25  A^. 


Fig.  14.  Transfer  function  ©  for  cold  chemistry  and  chemical  reaction  methods  displaying  the  sensitivity  of 
the  techniques. 
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Fig.  16.  Magnified  cold  chemistry  images  at  =  0.39  with  slow-side  seeding.  Flow  is  left-to-right  with 
the  high-speed  stream  on  top  and  the  imaged  region  is  Aa:  =4.8  cm  by  Ay  =5.0  cm.  Plotted  is  the  pure 
fluid  fraction  versus  distance  along  the  streamwise  and  transverse  lines  indicated. 


5.2.1  Diagnostic  Development 


Since  the  cold  chemistry  technique  used  to  quantitatively  measure  molecular 
mixing  is  relatively  new  (Paul  and  Clemens  (1993);  Clemens  and  Paul  (1995)), 
diagnostic  development  was  an  important  aspect  of  this  work.  In  this  method,  NO  is 
seeded  in  a  nitrogen  stream  that  mixes  with  an  unseeded  air  stream.  Fluorescence  from 
the  seeded  fluid  will  be  strongly  quenched  by  oxygen  for  fluid  that  has  been  mixed.  The 
quenching  occurs  by  molecular  collisions  in  a  time  scale  of  tens  of  nanoseconds 
compared  to  the  microsecond  fluid  mechanical  time  scale,  thus  mimicking  a  very  fast, 
fuel  lean  chemical  reaction  without  the  complication  of  heat  release.  Because  even  a 
small  amount  of  air  reduces  the  fluorescence  below  detection,  the  signal  is  approximately 
binary,  being  high  for  pure  fluid  and  low  for  mixed  fluid.  Therefore,  even  though  the 
probe  (pixel)  volume  may  be  significantly  larger  than  the  smallest  scales  of  mixing,  the 
cold  chemistry  method  returns  the  fraction  of  pure  seeded  fluid  present  in  the  probe 
volume.  Quantitative  statistical  mixing  quantities  can  be  derived  from  flip  experiments 
measuring  these  pure  fluid  profiles.  The  current  work  extends  the  technique  by 
combining  passive  scalar  and  cold  chemistry  methods  to  quantitatively  determine  the 
mixed  fluid  composition. 

Extension  for  Mixed  Fluid  Composition 

In  a  novel  extension  to  the  cold  chemistry  method,  the  mixed  fluid  composition 
can  be  determined  if  cold  chemistry  and  passive  scalar  measurements  are  combined. 
By  definition,  the  mean  high  speed  composition  is 

l(r\)  =  Pi(Ti)(l)  +  P2(ti)(0)  +  P„(Ti)^„  . 


where  |  is  the  mean  high-speed  mixture  fraction,  t|  is  the  transverse  coordinate,  Pj  is  the 
probability  of  pure  fluid  from  stream  i,  Pm  is  the  probability  of  mixed  fluid,  and  ^  is  the 
mean  composition  of  the  mixed  fluid.  Rearranging  gives  the  expression  for  the  mixed 
fluid  eomposition 


4m  (Tl)  = 


4(h)- Pi  (h) 

Pm(h) 


where  the  mean  composition  |  is  determined  by  the  passive  scalar  method  and  the  pure 
fluid  P]  and  mixed  fluid  profiles  Pm  are  determined  by  the  cold  chemistry  method. 
Essentially,  the  contribution  of  pure  fluid  to  the  mean  is  subtracted,  leaving  only  the 
contribution  of  mixed  fluid. 


21 


This  extension  to  the  cold  chemistry  measurement  is  significant  since  the  mixed 
fluid  composition  and  its  distribution  in  the  mixing  layer  is  critical  to  many  applications, 
such  as  supersonic  combustion. 

Signal  Modeling 

As  described  in  Island  (1997),  a  simple  two-level  model  was  developed  for 
weakly  and  strongly  quenched  NO  fluorescence  and  was  used  to  estimate  the  theoretical 
accuracy  of  the  cold  chemistry  method  and  compare  the  technique  to  chemistry-based 
mixing  diagnostics.  Figure  13  shows  the  complete  normalized  fluorescence  signal  for  the 
passive  scalar  and  cold  chemistry  methods  as  a  function  of  high-speed  mixture  fraction 
when  the  seeded  high-speed  stream  is  N2  and  the  unseeded  stream  is  air.  The  passive 
scalar  signal  decreases  linearly  with  mixture  fraction  due  simply  to  dilution.  The  cold 
chemistry  signal  decreases  non-linearly  due  to  quenching,  dominant  at  high  mixture 
fractions,  and  linearly  due  to  dilution,  dominant  at  low  mixture  fractions.  This  can  be 
seen  by  comparing  the  cold  chemistry  signal  with  the  curve  for  constant  mole  fraction  of 
NO  as  if  both  streams  had  been  seeded.  Because  of  the  strong  quenching,  the  signal  is 
reduced  ten-fold  for  6%  mixing,  twenty-fold  for  14%  mixing,  and  100-fold  for  50% 
mixing. 

Accuracy  and  Comparison  to  Chemical  Reaction  Techniques 

The  accuracy  of  the  cold  chemistry  method  can  be  examined  using  a  mixed  fluid 
transform  function.  From  the  PDF  of  mixture  fraction,  the  probability  of  mixed  fluid  at 
any  composition  is  defined  as 

where  the  arbitrarily  small  number  £  indicates  that  contributions  from  pure  fluid  from 
either  stream  are  excluded.  We  can  introduce  a  transfer  function  0  representing  the 
measurement  technique  such  that 

The  estimate  of  mixed  fluid  0n,  will  be  quite  good  for  transfer  functions  that  approach  the 
ideal:  zero  at  ^  =  0,  1  and  unity  elsewhere.  This  function  is  plotted  in  Figure  14  along 
with  the  transfer  function  for  the  chemical  reaction  method. 

The  cold  chemistry  technique  compares  favorably  with  the  chemical  reaction 
approach  for  practical  ({>.  One  measure  of  comparison  is  the  integral  of  the  transfer 
function  which  equals  0.93  for  cold  chemistry  and  0.80,  0.89,  and  0.94  for  ^  =  1/4,  1/8, 
and  1/16  chemistry.  Thus  as  a  guide,  the  cold  chemistry  sensitivity  is  equal  to  the 
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chemical  reaction  method  with  (j)  =  1/16.  The  absolute  error  depends  on  the  actual 
distribution  of  mixture  fractions,  i.e.  the  PDF  of  mixture  fraction.  Using  the 
incompressible  data  of  Frieler  (1992),  the  cold  chemistry  absolute  errors  are  estimated  to 
be  roughly  15%  at  the  edges  and  5%  near  the  center,  always  overestimating  pure  fluid. 
By  comparison,  the  chemistry  method  with  (|)  =  1/4  underestimates  the  amount  of  mixed 
fluid  by  40%,  7%,  and  9%  for  y/5  =  -0.4,  -0.2,  and  0,  respectively.  At  the  edge  of  the 
layer  the  significant  amount  of  mixed  fluid  near  the  freestream  composition  is  missed  at 
this  low  equivalence  ratio. 

5.2.2  Compressible  Mixing:  Imaging  Results 

PLIF  images  were  acquired  at  four  compressibilities  (Me  =  0.25,  0.39,  0.63,  and 
0.76)  and  three  Reynolds  numbers  for  each  compressibility  except  the  highest,  for  a  total 
of  10  cases.  Four  passive  scalar  images  are  shown  in  Figure  15  for  slow-side  seeding  at 
Me  =  0.39,  with  horizontal  and  vertical  lines  indicating  the  location  of  the  streamwise  and 
transverse  intensity  “cuts”  plotted  adjacent  to  each  image.  Since  these  are  quantitative 
images,  the  intensity  scale  represents  mixture  fraction  of  seeded  fluid.  Images  (a),  (b), 
and  (c)  are  clear  examples  of  Brown-Roshko  type  rollers.  As  seen  in  the  intensity 
profiles,  the  structures  are  often  “ramped”  in  ^  in  the  streamwise  direction  and  fairly 
uniform  in  the  transverse  direction.  This  is  evidence  of  entrainment  occurring  by  2D 
motions  with  the  structure  core  being  biased  in  mixture  fraction  to  the  side  from  which  it 
entrains  pure  fluid.  For  the  images  shown,  the  mixture  fraction  in  the  roller,  core  is  about 
1/2,  but  structure-to  structure  variation  and  streamwise  ramping  result  in  a  smooth  time- 
average. 

Figure  16  shows  typical  cold  chemistry  images  for  the  same  low  compressibility 
case.  The  seeded  low-speed  stream  is  nitrogen  and  the  high-speed  stream  is  air  so  that 
significant  fluorescence  quenching  occurs  in  the  mixed  fluid  by  oxygen  molecules.  The 
cold  chemistry  images  are  nearly  binary,  with  the  pure  seeded  fluid  white  and  mixed  fluid 
black,  in  contrast  to  the  many  gray  levels  in  the  passive  scalar  images  of  Figure  15.  In 
general,  the  2D  structure  of  the  layer  is  harder  to  distinguish  from  the  pure  fluid  images. 
As  expected,  however,  many  of  the  images  (e.g.,  (b)  and  (c))  strongly  suggest  that  pure 
fluid  from  stream  2  is  entrained  on  the  upstream  edge  of  the  2D  roller  (the  tongues  of 
pure  fluid)  and  is  well  mixed  in  the  roller  core.  As  seen  in  the  intensity  cuts,  transverse 
profiles  though  the  mixing  layer  in  regions  not  associated  with  presumed  2D  structure, 
such  as  Figure  16  (d),  show  a  rapid  fall-off  in  the  pure  fluid  with  transverse  distance. 
There  is  a  relatively  sharp  interface  between  pure  and  mixed  fluid  in  these  regions  and 
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the  mixing  layer  contains  very  little  pure  fluid.  However,  even  in  these  regions,  a  thin 
band  of  low  signal  (about  10%)  along  the  pure  fluid  interface  can  just  be  discerned  in  the 
images  and  clearly  seen  in  the  transverse  intensity  profiles.  This  “halo”  is  most  likely 
pure  fluid  present  at  scales  below  the  resolution  of  the  imaging  system.  From  the  halos, 
we  assume  that  pure  and  mixed  fluid  is  being  stirred  by  the  action  of  small-scale 
turbulence  along  the  interface. 

The  following  conclusions  regarding  the  effects  of  compressibility  may  be 
generalized  from  the  entire  image  set: 

1)  As  seen  by  previous  researchers,  the  passive  scalar  visualizations  show  a  change  in 
large-scale  structure  with  increased  compressibility.  While  somewhat  variable,  the  low 
compressibility  image  sets  are  dominated  by  two-dimensional  rollers  exhibiting 
streamwiseramping  and  cross-stream  uniformity  in  mixture  fraction.  At  the  higher 
compressibilities,  the  large-scale  structure  is  highly  variable  with  very  few  images 
exhibiting  strongly  two-dimensional  structure. 

2)  Similarly,  the  cold  chemistry  visualizations  of  pure  fluid  indicate  a  change  in  structure 
with  compressibility.  At  low  compressibility,  “tongues”  of  pure  fluid  appear  to 
correspond  with  entrainment  of  high-  and  low-speed  pure  fluid  along  the  leading-  and 
trailing-edges,  respectively,  of  two-dimensional  structures  and  extend  deep  into  the 
mixing  layer.  With  higher  compressibility,  the  pure  fluid  intrusions  are  less  periodic  and 
have  a  reduced  transverse  extent,  indicating  that  entrainment  is  less  two-dimensional. 

3)  The  cold  chemistry  images  frequently  show  a  “halo”  of  low-signal  signal  level 
adjacent  to  pure  fluid  regions.  In  the  low  compressibility  case,  such  regions  extend  well 
into  the  mixing  layer  and  can  be  found  in  the  (presumed)  core  regions  of  two-dimensional 
structures.  At  higher  compressibility,  the  halo  regions  are,  in  general,  less  prevalent  and 
more  confined  to  the  mixing  layer  edge. 

5.2.3  Compressible  Mixing:  Statistical  Results 
Mean  CompositiQn 

The  profiles,  of  mean  composition  or  high-speed  mixture  fraction  ^  vs.  the 
normalized  transverse  coordinate,  r\,  are  shown  in  Figure  17.  The  Me  =  0.63  and  0.79 
passive  scalar  results  are  corrected  for  air  contamination  and  should  be  interpreted 
cautiously  due  to  the  uncertainty  associated  with  the  correction.  The  profiles  of  Figure 
17  show  fairly  strong  collapse  to  a  single-inflection  point  profile.  No  clear  trend  with  Mg 
or  Re  is  evident  within  the  uncertainty  of  the  measurement,  and  we  conclude  that  ^(ri)  is 
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not  strongly  influenced  by  compressibility  in  the  range  studied.  In  the  present  work,  the 
mixing  layer  scalar  thickness  5  is  defined  as  the  transverse  distance  between  the  5  and 
95%  points  of  the  mean  composition  profiles  In  agreement  with  numerous  studies, 
there  is  a  strong  decrease  in  mixing  layer  thickness  normalized  by  its  incompressible 
value  at  the  same  density  and  velocity  ratio  with  increased  compressibility. 

Pure  Fluid 

The  probability  of  pure  seeded-side  fluid  Pi(r|)  and  P2(r|)  is  plotted  in  Figure  18. 
This  quantity  is  the  normalized  mean  fluorescence  signal  of  the  cold  chemistry  flip 
experiments.  The  measurements  indicate  that  the  pure  fluid  transverse  profiles,  when 
scaled  on  rice  based  on  the  cold  chemistry  measurements  themselves,  are  independent  of 
Reynolds  number  and  compressibility,  strongly  collapsing  to  a  single-inflection  point 
curve. 

The  pure  fluid  profiles  are,  however,  shifted  and  stretched  with  respect  to  the 
mean  composition  (passive  scalar)  midpoint  and  thickness.  To  demonstrate.  Figure  19 
displays  the  probability  of  pure  fluid  versus  t]  based  on  the  corrected  passive  scalar 
profiles.  In  Figure  19(a),  the  probability  of  high-speed  fluid  at  a  given  transverse  location 
is  seen  to  decrease  strongly  with  compressibility.  No  such  effect  is  indicated  in  Figure 
19(b)  along  the  low-speed  interface  for  the  low-speed  fluid.  Thus,  we  expect  that  the 
mean  mixed  fluid  composition  will  be  somewhat  higher  along  the  high-speed  interface 
and  relatively  unchanged  along  the  low-speed  interface  with  increased  compressibility. 

Since  instantaneous  images  were  acquired,  the  root  mean  square  (RMS) 
fluctuations  in  the  pure  fluid  profiles  Pi(r|)  can  be  calculated.  Figure  20  shows  two 
interesting  trends  with  increased  compressibility:  1)  the  peak  fluctuations  shift  towards 
the  edges  of  the  mixing  layer  for  both  P]  and  P2,  and  2)  the  fluctuation  peak  narrows, 
mostly  on  the  high-speed  edge.  Additionally,  the  fluctuation  peak  is  generally  broader 
along  the  low-speed  edge  than  along  the  high-speed  interface. 

Probability  of  Mixed  Fluid 

From  the  pure  fluid  profiles,  the  probability  of  mixed  fluid  profiles  Pm(il)  are 
generated  and  plotted  in  Figure  21.  The  peak  value  of  Pm  is  always  less  than  unity  even 
at  the  midpoint,  in  agreement  with  other  researchers  using  chemically  reacting  and  cold 
chemistry  methods  (Hall  (1991);  Clemens  and  Paul  (1995)),  although  at  the  higher 
compressibilities,  the  chance  of  finding  mixed  fluid  is  better  than  95%  for  a  fairly  broad 
region  of  the  mixing  layer.  With  increasing  Me  and  Re  the  chance  of  finding  mixed  fluid 
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Fig.  17.  /?e-averaged  mean  composition  profiles:  fast-  and  slow-side  seeding  combined  and  corrected  for  air 
contamination. 
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RMS  of  Pi 


Fig.  21.  Probability  of  mixed  fluid  profiles  from  the  cold  chemistry  method,  with  T|  scaled  on  the  conected 
passive  scalar  mean  composition. 
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Fig.  22.  Mixing  efficiency  as  a  function  of  (a)  and  (b)  Re  for  the  present  and  other  resolution-insensitive 


measurements. 


for  air  contamination. 
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becomes  biased  towards  the  high-speed  interface  suggesting  that  the  mixed  mean 
composition  near  the  high-speed  edge  will  increase  with  compressibility. 

The  mixing  efficiency  dm/d  denotes  the  fraction  of  fluid  in  the  layer  which  is 
molecularly  mixed  and  is  defined  as  the  integral  of  the  probability  of  mixed  fluid  profile 
Pm(Tl).  The  mixing  efficiencies  determined  from  the  10  cold  chemistry  flip  experiments 
are  plotted  versus  convective  Mach  number  in  Figure  22(a)  and  Reynolds  number  in 
Figure  22(b).  As  displayed  in  Figure  22(a),  the.re  is  an  apparent  trend  of  moderately 
increasing  mixing  efficiency  with  increasing  compressibility,  with  6m/5  =  0.50,  0.55, 
0.60,  and  0.62  at  Me  =  0.25,  0.39,  0.63,  and  0.76,  respectively.  Hall  (1991),  using  a 
chemical  reaction  method  in  planar  mixing  layers,  found  mixing  efficiency  to  decline  by 
15%  between  Me  =  0.5  and  0.9;  however,  his  technique  suffers  from  relatively  poor 
equivalence  ratio  («t)=l/4)  and  potential  finite  rate  chemistry  effects  due  to  short  fluid 
mechanical  time  scales  in  the  compressible  mixing  layer. 

The  apparent  trend  of  increased  mixing  efficiency  with  compressibility  seen  in  the 
present  experiments,  however,  may  be  attributable  to  Reynolds  number  effects.  As  seen 
in  Figure  22(b),  the  present  results  show  a  strong  correlation  of  increased  mixing 
efficiency  with  Re,  approximately  20%  per  decade.  Except  for  a  single  point  at  low 
compressibility,  each  of  the  experiments  conducted  at  a  given  compressibility  shows  an 
increase  in  5m/5  with  Re.  Additionally,  the  highest  compressibility  point  has  nearly  the 
same  Re  as  one  of  the  Me  =  0.63  cases  and  has  a  similar  mixing  efficiency,  despite  a  21  % 
increase  in  Me.  The  strong  correlation  between  and  Re  seems  to  largely  account  for 
the  change  in  efficiency  seen  with  compressibility,  and  we  therefore  conclude  that  mixing 
efficiency  is  essentially  independent  of  compressibility. 

Mixed  Mean  Composition 

As  discussed  earlier,  the  mixed  mean  composition  i.e.,  the  average  high-speed 
mixture  fraction  of  fluid  which  has  been  mixed,  is  calculated  by  subtracting  the 
component  of  pure  fluid  as  determined  by  the  cold  chemistry  measurement  from  the 
mean  composition  determined  by  the  passive  scalar  technique.  Figure  23  displays  the 
results.  Increased  compressibility  appears  to  shift  the  mixed  mean  away  from  a  preferred 
value  and  toward  a  more  continuous  or  freestream-biased  profile.  This  trend  is  associated 
with  compressibility  and  not  density  ratio  because  the  highest  and  lowest  compressibility 
cases  have  approximately  the  same  density  ratio.  Note  that  most  of  the  change  appears  to 
occur  along  the  high-speed  interface  and  that  even  under  compressible  conditions  little  of 
the  mixed  fluid,  on  average,  is  comprised  of  more  than  80%  freestream  fluid.  Again, 
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these  results  should  be  regarded  cautiously  since  the  passive  scalar  component  is  subject 
to  air  contamination  in  the  more  compressible  cases.  The  present  results  generally 
support  the  reacting  results  of  Miller  et  al.  (1997)  and  would  explain  more  intense  OH  in 
compressible  cases. 

5.2.4  Mixing  Enhancement 

In  the  mixing  enhancement  study,  sub-boundary  layer  disturbances  were 
introduced  by  applying  small  pieces  of  Scotch  tape  having  a  thickness  of  60  microns  to 
the  high-speed  (supersonic)  side  of  the  splitter  tip  as  shown  in  Figure  24.  The  geometric 
parameters  varied  were  1)  shape  of  the  disturbance  (from  sharp  triangles  to  squares)  as 
determined  by  (j),  2)  number  N  of  disturbances,  3)  obliquity  0  of  an  array  of  disturbances, 
4)  offset  Lo  from  the  splitter  tip,  and  5)  thickness  t  (height  from  the  splitter  tip  surface)  of 
the  disturbance.  The  thickness  was  varied  by  layering  the  Scotch  tape  and  ranges  from  5 
to  50%  of  the  calculated  splitter  tip  boundary  layer  displaeement  thickness  899. 
Additionally,  two-dimensional  disturbances  with  single  or  double  spanwise  trips,  as 
shown  in  Figure  24(d),  were  tested  for  smooth  and  rough  surfaces.  For  each  geometry, 
time  resolved  end,  side,  and  plan  view  visualizations  of  mixed  fluid  were  obtained  with 
the  planar  laser  Mie  scattering  (PLMS)  technique.  The  product  formation  method  was 
used  whereby  alcohol  droplets  condense  in  mixed  fluid  and  high  scattering  signal 
consequently  marks  mixed  fluid. 

5.2.5  Mixing  Enhancement  Results:  Images 
Schlieren 

For  reference,  Toepler  schlieren  side  view  images  are  shown  in  Figure  25  for  the 
Me  =  0.76  case  with  and  without  a  disturbance.  The  disturbance  geometry  consists  of 
five  discrete  equilateral  triangles  located  on  the  end  of  the  splitter  tip,  with  a  thickness 
one-half  of  the  calculated  boundary  layer  displacement  thickness. 

For  the  unperturbed  case.  Figure  25(a)  shows  the  shock  structure  and  mixing  layer 
for  the  first  12  cm  of  the  test  section.  The  shock  originating  upstream  of  the  splitter  tip  is 
a  result  of  the  (very)  slight  mismateh  between  the  nozzle  wall  and  upper  guidewall,  and  a 
second  shock  is  generated  from  the  flow  expansion  at  the  end  of  the  splitter.  The  mixing 
layer  grows  approximately  linearly  from  a  virtual  origin  upstream  of  the  splitter  tip.  In 
the  elose-up  shown  in  Figure  25(b),  the  high-speed  boundary  layer  can  just  be  discerned, 
especially  in  the  average  view. 
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For  the  perturbed  case.  Figure  25(c)  shows  a  dramatic  increase  in  the  mixing  layer 
thickness  and  increased  growth  rate,  with  5  about  twice  that  of  the  unperturbed  case. 
Since  Toepler  schlieren  integrates  across  the  mixing  layer,  it  cannot,  however,  distinguish 
between  the  entire  layer  thickening  or  the  presence  of  spanwise  convolutions.  Most  of 
the  increased  growth  occurs  along  the  low-speed  interface  and  the  virtual  origin  shifts 
towards  the  splitter  tip.  Comparing  to  Figure  25(a),  a  new  shock  is  evident  that  originates 
upstream  of  the  splitter  tip  and  penetrates  the  mixing  layer  after  reflection  from  the  upper 
wall.  From  the  close-up  average  in  Figure  25(d),  this  shock  can  be  seen  to  originate  from 
the  upstream  edge  of  the  disturbance.  Although  this  is  the  thickest  disturbance 
investigated.  Figure  25(d)  shows  it  well  within  the  (thermal)  boundary  layer. 

PLUS 

For  brevity,  planar  laser  Mie  scattering  images  at  the  Me  =  0.63  flow  condition  are 
presented  here  (see  Figure  26)  only  for  the  same  geometry  as  the  Schlieren  images.  Flow 
is  left  to  right  with  instantaneous  and  average  (a)  plan  views,  (b)  end  views,  and  (c)  side 
views.  The  plan  views  of  Figure  26(a)  clearly  indicate  streamwise  structure  which  is 
well-localized,  as  seen  in  the  average.  Even  with  the  thin  disturbance  (t/599  =  5%)  used  in 
the  end  view  imaging,  significant  spanwise  convolutions  can  be  discerned.  The  bottom 
end  view  image  is  included  as  a  typical  example  that  suggests  the  presence  of  counter¬ 
rotating  vortex  pairs  associated  with  the  discrete  disturbances.  To  be  objective,  there  is 
considerable  variation  in  the  end  views  in  all  cases,  but  images  suggesting  counter¬ 
rotating  vortices  are  both  more  frequent  and  more  consistent  for  discrete  3D  disturbances 
compared  to  the  unperturbed  case.  The  side  views  of  Figure  26(b)  also  indicate  a  thicker 
mixing  layer  for  this  case  relative  to  the  unperturbed  case. 

5.2.6  Mixing  Enhancement  Results:  Statistics 

The  mixing  layer  thickness  was  measured  directly  from  the  side  and  end  view 
PLMS  images  to  examine  the  effect  of  different  disturbance  geometries.  The  most 
significant  results  are  described  below  and  further  detail  can  be  found  in  Island  (1997). 

Shape 

The  effects  of  disturbance  shape  on  mixing  layer  thickness  is  shown  in  Figure  27 
which  plots  the  normalized  thickness  as  a  function  of  shape  parameter  Both  the  end 
view  measurements  (top)  and  side  view  measurements  (bottom)  suggest  that  the  optimal 
shape  is  an  equilateral  triangle  (j)  =  30  degrees.  The  most  conclusive  data  comes  from  the 
thicker  disturbances  in  the  side  view  measurement  where  the  thickness  is  20%  greater  for 
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Fig.  25.  Schlieren  images  at  Me  =  0.76  for  (a)  unperturbed  wide  view,  (b)  unperturbed  close-up,  (c) 
perturbedwide  view,  and  (d)perturbedclose-up.  Shown  are  (left)  instantaneous  and  (right)  averagedimages. 


Fig.  26.  Typical  instantaneous  and  average  PLMS  images  for  (a)  plan,  (b)  side,  and  (c)  end  views  for  the 
case  3  geometry.  The  disturbance  thickness  is  =  25%  for  the  plan  and  side  views  and  =  5%  for  the 
end  views. 


(j)  (degrees) 


Fig.  27.  Normalized  mixing  layer  thickness  5/5o  with  different  shape  elements  measured  from  (top)  end 
views  and  (bottom)  side  views.  As  shown,  the  shape  is  specified  by  the  angle  relative  to  the  streamwise 
direction,  so  that  a  square  is  90°  and  an  equilateral  triangle  is  30°. 


Fig,  28,  Normalized  mixing  layer  thickness  6/5o  at  differentoffsets  Lq/L  from  the  splitter  tip,  as  defined  in 
the  right  figure.  Determined  from  side  and  end  views. 
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Fig.  29.  Normalized  mixing  layer  thickness  5/5o  for  2D  and  3D  triangular  perturbations  of  different 
thickness.  The  disturbance  thickness  5/899  is  normalized  on  the  high-speed  boundary  layer  displacement 
thickness.  Determined  from  side  and  end  views. 
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an  equilateral  triangle  compared  to  the  more  and  the  less  acute  triangles  tested.  Note  that 
the  flow  area  blockage  is  the  same  for  all  the  shapes  tested  so  we  expect  roughly  the  same 
pressure  losses  associated  with  each  shape.  Square  shapes  ((|)  =  90  degrees)  were  found 
to  be  the  least  effective  at  enhancing  thickness  even  though  the  images  do  show  a 
moderate  increase  in  streamwise  structure. 

Offset 

The  impact  of  offsetting  the  disturbances  from  the  splitter  tip  is  shown  in  Figure 
28.  As  also  suggested  by  the  images,  the  largest  improvement  in  6  occurs  for 
disturbances  placed  directly  at  the  splitter  tip.  Again,  the  most  dramatic  evidence  comes 
from  the  thicker  disturbances  measured  in  side  view  which  show  a  45%  improvement  in 
5  for  disturbances  with  zero  offset  as  opposed  to  almost  no  improvement  for  an  offset  of 
Lo/L  =  1/2. 

Disturbance  Thickness  and  2D  V5.  3D  Effects 

Figure  29  plots  the  improvement  in  layer  thickness  versus  the  disturbance 
thickness  for  discrete  triangular  and  2D  disturbances.  The  triangular  disturbances 
dramatically  increase  the  layer  thickness  with  a  30  and  60%  thickening  for  a  disturbance 
thickness  of  15  and  50%  of  the  boundary  layer,  respectively.  Equally  notable  is  the 
negative  result  for  the  2D  and  double  2D  spanwise  disturbance  which  actually  slightly 
decrease  the  layer  thickness,  despite  having  equivalent  flow  blockage  as  the  triangles. 
Smooth  (Scotch  tape)  and  rough  (400  and  120  grit  sandpaper)  2D  disturbances  show 
similar  results. 

5.2.7  Mixing  Enhancement  Results:  Quantitative  Cold  Chemistry  Measurements 

Cold  chemistry  measurements  were  made  for  the  Me  =  0.63  flowfield  for  the  most 
promising  enhancement  geometry  (five  equilateral  triangles  at  the  splitter  tip)  having 
thickness  t/599  =  25%.  The  principal  findings  are:  1)  the  enhanced  case  shows  a  48% 
increase  in  mixing  layer  thickness  compared  to  the  cold  chemistry  thickness  of  the 
unperturbed  case  at  this  compressibility  and  Reynolds  number;  2)  the  probability  of  pure 
fluid  curves  are  quite  similar  between  the  perturbed  and  unperturbed  cases.  The  shape  of 
these  curves  now  appears  quite  robust  with  near  perfect  collapse  for  a  ten-fold  increase  in 
Re,  significant  changes  in  Me,  and  unperturbed  and  moderately  perturbed  mixing  layers; 
3)  the  pure  fluid  probabilities  show  greater  RMS  fluctuation  along  both  interfaces  for  the 
enhanced  case;  4)  for  the  perturbed  case,  the  probability  of  finding  mixed  fluid  is 
essentially  unity  in  the  middle  10%  of  the  mixing  layer,  compared  to  about  a  95% 


41 


F 


probability  in  the  unperturbed  Me  =  0.65  case.  Even  more  than  the  unperturbed 
compressible  cases,  Pm  is  biased  towards  the  high-speed  side.  In  contrast,  along  the  low- 
speed  edge  there  is  slightly  less  chance  of  finding  mixed  fluid  in  the  enhanced  layer;  5) 
the  measured  mixing  efficiency  is  5m/5  =  0.64  compared  to  0.60  for  the  unperturbed  layer 
at  the  same  Me  and  Re,  a  7%  increase  which  is  within  the  uncertainty  of  the  cold 
chemistry  measurements.  The  enhanced  case  has,  therefore,  approximately  57%  more 
mixed  fluid,  mostly  because  of  the  thicker  layer. 

The  quantitative  mixing  enhancement  results  are  particularly  significant  because 
they  show  that  a  qualitatively  thicker  layer  can  have  a  correspondingly  greater  amount  of 
molecularly-mixed  fluid.  That  is,  the  mixing  efficiency  of  the  layer  does  not  decrease  as 
a  result  of  the  thickness  enhancement  and  greater  bulk  entrainment  thereby  implies 
greater  entrainment  into  the  mixed  state,  which  is  the  relevant  issue  for  most  applications. 
This  result  is  not  necessarily  true  and  does  not,  of  course,  universally  apply  to  other 
enhancement  strategies  which  might  thicken  the  mixing  layer  by  increasing  the  amount  of 
pure  (unmixed)  fluid  in  the  layer.  However,  the  constancy  of  mixing  efficiency  in  this 
case  bodes  well  for  the  applicability  of  the  large  body  of  qualitative  mixing  enhancement 
results  in  the  literature. 
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